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Abstract 

We use the Peierls-Yoccoz projection method to study the motion of a relativistic system of 
nucleons interacting with sigma and omega mesons, generalizing a method developed for the alpha 
particle. The nuclear system is described in a mean-field Hartree approach, including explicitly the 
meson contribution. The formalism is applied to ^He, ^^O and ^"^Ca. The center-of-mass correction 
makes the system too much bounded. It turns out that a new set of model parameters is needed 
when the center-of-mass motion is consistently treated with respect to the traditional approaches. 
An appropriate refitting of the model brings the radii and binding energies to reasonable values 
for the oxygen and calcium. 
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PACS number(s): 21.10.-k, 21.10.Dr, 21.60.-n 



I. INTRODUCTION 



Relativistic models for finite nuclei, with nucleons and mesons are usually treated in the 
Hartree or Hartree-Fock approximations. In these approximations the total linear momen- 
tum is no longer a conserved quantity and the spurious center of mass (CM) motion gives 
rise to unphysical contributions in the calculated nuclear observables. As already discussed 
for non-relativistic nuclear models, using the same kind of mean-field approximation, the 
correct treatment for the CM motion introduces a modest modification in the total binding 



energy for intermediate mass nuclei, but relatively large contributions in other observables, 
e.g., charge distributions and spectral functions In relativistic treatments, the CM 



correction is up to now limited to the harmonic approximation for the energy or to the sub- 



traction of „A. from the total energy, where Pa is the total nucleonic momentum operator. 



2AM 



M is the nucleon mass and the mean-value is taken using the Hartree self-consistent state 



for A nucleons. 

More recently 3|, the CM energy correction was estimated within the a — ui model. 



using the Peierls-Yoccoz projection procedure 



J], for N = Z spherical nuclei within the 



Hartree approximation. Although the correction obtained in this way is of the same order 



of magnitude of the harmonic approximation, only the nucleonic degrees of freedom were 
taken into account in that calculation. In reference 5|, a formalism has been developed to 



include the mesonic degrees of freedom in the CM projection within a — u models and an 
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application was then made to the ^He nucleus. In the present paper we generalize the results 
obtained in 51] in order to extend the calculations to heavier spherical nuclei, allowing us 

;o draw more systematic conclusions. In section II, we review the main results from Ref. 

5|. Then, in section III, the linear momentum projection within the model is presented and 
the total energy functional is worked out. Since most of the model parametrizations within 
a — uj models are based on fits to the experimental data of both binding energy and charge 
radius, we take the same point of view. The nuclear charge mean square radius is discussed 
in section IV. The numerical results for ^He, ^^O and ^°Ca are shown and discussed in section 
V. Finally the conclusions and perspectives are summarized in section VI. 

II. THE MEAN-FIELD HAMILTONIAN WITH MESONS 
AS COHERENT STATES 

In this section we summarize the main aspects of the relativistic nucleon-meson models of 
nuclei. The model used in this work is restricted to a and uj mesons, without self-interactions 
(the inclusion of such interactions is straightforward using the method presented in this 
paper). 

The Lagrangian density for a system of nucleons interacting with sigma and omega mesons 



reads 



T T free i j free i r free i r int , r int / -i '\ 
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where N denotes the nucleon and cr, oj the mesons. The Lagrangians for the free fields are: 



T free 



--[mla\x)-d,aix)d^a{x)], 
-^F^,{x)F^^{x) + ^mlu;,{x)u;''{x), 



(2) 
(3) 
(4) 



where 



Fiuu = d^uj^{x) - d^uj^{x) , 



M is the rest mass of the nucleon, and rria and rriu, are the meson masses. The sigma and 
omega fields are denoted respectively by a{x) and u;'^{x), and the nucleon field by tlj{x). The 
interaction parts of the Lagrangian are 



^nL = gai^{x)cr{x)ijj{x) 



(5) 
(6) 



Prom the above Lagrangian density one derives the following Hamiltonian density: 



(7) 



where the fermionic term is 



9u> ^ ^ 



Hn = r(x) \-ia-V + (3[M - g^a(x)] - g^a ■ u{x) + • \ ^P{x) 



mf 



(8) 



and the free meson terms are: 



P ■ P + 



[V-P^ 



+ (V X ujf + m^tJ^ 



(9) 



5|, but the procedure can be 



n„ = ]^[Pl + Va-Va + mla^] . (10) 

Note that we have used the definitions = with i = 1,2,3 and = doc. The 
quantization of the model follows the usual procedure described, for instance, in [7|. The 
most important steps of the quantization are described in 
be outlined here by saying that, for the canonical quantization of massive vector fields, one 
cannot use the field uq, because its canonical conjugate field is identically zero. For conserved 
four- vector sources (as it is the case of the nucleon vector current) the four-divergence of a;^ 
is zero, and therefore one can use the full Klein-Gordon equation for ujq to write this field 
in terms of the divergence of the conjugate field of the spatial components tu* and of the 
zero-th component of the vector current, gcu'ip'^'ip (see eq. (12) in Ref. [5|). The Hamiltonian 
is then built in the usual way by using only the spatial components, u!\ and their respective 
conjugate fields, P^. The two-body contact term in ([8]) arises from the quadratic term in uq 
in the Lagrangian. 

The nucleon field operators can be expanded as 

^(x) = J2 ^"(^1 e-^^'^* b^ + Y, e'^"* 4 (11) 

a a 

tlj^x) = J2 4(^1 e'"""* bi + Y, ^1(0 e-'^'^* do, , (12) 

a a 

where Ua{r) and Va{r) form a complete set of Dirac spinors in the coordinate space, and ba 
and are the creation and annihilation operators of a nucleon in the state a. By da and 
dl^ we denote the creation and the annihilation operators for the anti-nucleons in the state 
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a. Similarly, the a meson field may also be expanded in the following form: 



c{k) e'^-^>c^(A;) e-''^'^^ 



(13) 



The omega field expansion, considering longitudinal and transverse waves relative to the 



wave vector /c, reads: 



kai{k) + ^ et{k)at{k) 



t=l,2 



e^^^-^'+h.c \ . 



(14) 



All creation and annihilation operators {c^c^)^{ai^a\) and (at, a|) obey canonical boson 



commutation relations, and we have introduced the frequencies Ua- — + k'^ and 



\/vn^ + k'^. Using the above expansions, the free meson field Hamiltonians can be 



cast in the form: 



E^^ Jd^rHa^ J d^kuja{k)c'{k)c{k), 



(15) 



and 



= / d^ku^k) 



a}{k)ai{k) + J2 4{k)at{k) 

t=l,2 



(16) 



The nucleus state is assumed to be described by |V' >= \A > \a > \u; >, with 1^4 > 
representing an A fermion Slater determinant with the lowest energy states occupied (valence 
or no-sea approximation), i.e.. 



A>^blbl---bi^ |0>, 



(17) 



where cti, . . . , are sets of single-particle quantum numbers and | > is the bare vacuum. 
As it is usual in o" — a; models, we work in the valence approximation, which means that the 



polarization of the negative energy single particle states is neglected. This is also a common 



approximation in quark-meson chiral soliton models, such as the hnear sigma model or 
the chromodielectric model 8|. The approximation is even more justifiable here since the 



binding energy per nucleon is small compared with the rest mass of the nucleon. However, 
this might not be the case in other approximations (see, for instance, Ref. 9|). 



In the above product state, \a > represents a coherent state describing the a mesons and 
\uj > a coherent state describing the u mesons. For instance, for the a meson cloud: 



|cr >= Ncr exp 



10 > 



(18) 



with c{k)\a >= r]{k)\a > and, from the normalization of the state 



exp 



d^k\T]{k)\'' 



(19) 



We now enforce the mean value of the a field operator in the coherent state to be equal to 
the potential obtained in the mean-field Hartree approximation, i.e., we demand 



< a\a\cr >= (po{r) 



(20) 



for a spherical symmetric potential. This condition allows us to determine, in an unique 
way, the function r]{k) in f|T8|) . Exploiting the spherical symmetry of the scalar potential, 
one finds: 



TT 



(21) 



where jo is the spherical Bessel function of zeroth order. A similar procedure can be carried 
out for the co meson field. In this case. 



\uj >= exp 



d^k [Qi{k)a]{k) + J2 ^t{k)al{k) 



10 > 



(22) 



and, using again the normalization and the properties of the vector potential in the Hartree 

Ml 

dr 



approximation, such as < cuplu; >= 0, < culcu^lcu >= uJo{r) and < cul-P^^lct; >= f — 



finds: 



n,(k) = [ .rr^Mtr)'-^, (23) 

and ^t{k) = 0. If we now take the states defined in flTTI) . ffTS]) and (1221) and calculate the 
mean value of the Hamiltonian obtained from ([7]), we exactly recover the nucleus energy 
obtained in the usual Hartree approximation. 

Let us stress that the coherent state is a multiparticle state and the description of meson 
clouds by coherent states introduces many body correlations. 

III. THE CENTER-OF-MASS APPROXIMATE PROJECTION 

Next, we want to obtain the center-of-mass (CM) correction to the energy using the model 
described in section II. It is well known, from the nuclear many-body theory, that mean-field 
approximations break translational invariance (see Ref. ^) and that the broken symmetry 
can be recovered by applying the Peierls-Yoccoz projection to the symmetry-breaking state. 
The projection operator 

Pp-= J exp[i{P - p) ■ a\d^d , (24) 

exhibits the property 

Vp-Vp-, = 6{p-p')V^. (25) 

In fl24|) . P = + P„ + is the total linear momentum operator and p the correspond- 
ing eigenvalue. Our approach consists in assuming that the model state representing the 



physical nucleus is obtained by projecting the product mean-field Hartree state onto a zero 
momentum {p = 0) state (the procedure is known as projection after variation). Since the 
Hamiltonian H = J d^rTi, with 7i given by ([7]), commutes with the projection operator, we 
may write the total energy, already corrected for the CM spurious motion, as 



E,=o = , ^ ' " • (26) 







ip > 


< i) 


1 "P^ - 1 

1 ' p=0 1 





We emphasize that, in the valence approximation, the projection operator acts on the mesons 
and on the positive energy fermions. The vacuum single-particle states are unperturbed and 
the vacuum is invariant under translations, so that the shifted states have the same energy 
as the unshifted ones in the so-called variation before projection method 10| that we are 
using here. 

In order to compute the projected energy let us first consider the norm overlap: 

< ^ I ^p-=o I ^ >= J da< a|e^^''-'^|a >< Lu\e'^--'^\uj >< A|e'^^-'^|v4 > , (27) 
and begin with the a field contribution. Its norm overlap reads: 

Ar^(a) =<a|e^^'^-"|(T>=exp|47ry" dk k^\r]{k)\^[jo{ka) - 1]^ , (28) 
where ri{k) is defined by equation (12T|) . We then find: 

N^ia) = exp |47r J dkk' ^ '^o(fe)bo(fca) - 1] j , (29) 

where 

4)o{k) = [ drr^3o{kr)(t)o{r). (30) 
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Similarly, for the uj meson norm overlap 

r 47r /■ irr? + Fl^/^ ^ 

= exp |— y ^ " ^ ^ cDg(fc) [jo(A;a) - 1] | , (31) 

where 

^o{k) = j drr'^ jo{kr)uJo{r). (32) 

The calculation of the fermionic part of the norm overlap is more involved, and we just quote 
here the main result in a compact form: 

NA{a) =< Ale'^^-'^IA >= det B , (33) 

where the B matrix is defined by 

B^p =< a\f3{a) > . (34) 

Each label (a and stands for the set of particle quantum numbers {n, m) as well as 
for the isospin projection quantum number necessary to classify the state. The ket |/3(a) > 
means a single-particle (four-component) state for which the spatial coordinate r is changed 
to r + a. 

Next, we move our attention to the energy kernel calculation. The total Hamiltonian is 
written in the form: 

H = Hm + H„ + H^, (35) 

where the first term contains the free fermion part as well as their interaction with the a — uu 
mesons. The second and third terms are given by equations ( fT5l) and ( fT6|) and represent the 
free mesonic terms. Let us consider the free a field energy kernel. Using equations ( fT8|) . ( l2T|l 
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and the result: 



cr(a) >= e "'"Icr >= iV<^exp 



dkr]'{k)b{k) 



|0>, 



(36) 



with 7] (k) = ri{k) e , we obtain: 



Eaia) =< a\HJa(a) >- 



dkk\mi + e)<i)i{k)]Q{ka) 



(37) 



For the free uj meson energy kernel, a similar analysis leads us to the following result: 



eJa) =< uj\Hjuj(a) >- 



2m2 



dkk\mi + k^)u'Q{k)jo{ka) 



NM ■ (38) 



For the fermionic part of the energy kernel, it is more convenient to rewrite the corresponding 
original Hamiltonian. From equation ([8]), in the Hartree mean-field, we can read off the 
fermionic Hamiltonian written in second quantized form: 



(39) 



with 



"■a/3 



j drulif) |-ia . V + /3 [M - g^a{x)] + ■ i^j M^(r) 



(40) 



and 




dfdf 'ul^{r)ut{f ')-^6{f— f ')u^{r)us{r '). 



m 



(41) 



In the above equations, Ua,i3 represents the Dirac single-particle spinor, which we choose 
to be the Hartree mean-field solution. Observing now that, the field should obey the 



Klein-Gordon equation: 



(42) 
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and that V ■ -FI, = —V'^ujQ(r), we may rewrite the one-body part in fl39l) as: 

/i('^ = W-^PB(r), (43) 
ml 

with: 

huFAUa = (^aUa- (44) 

We are now in position to perform the calculation of the fermionic part of the energy 



kernel, which reads 

eN{a) =< A\HNe'^^-'^\A >=J2^^ < A\V^'^e'^^-'^\A > + < A\h^' 2)giP^-a|^ ^ ^ ^45) 

a 

where we have defined V^^^ = -^pB{r)- The second and third terms in equation (l45l) can 



then be obtained with the help of the well-known results (see, e.g., Ref. 



< A I V^^^e'^^-'^ I A >= iV^(a) J]] < « I ^^^^ I /5(«) > > (46) 

Q/3 



and 



where the exchange term has been neglected. Putting everything together, we finally obtain 
the total nucleus energy corrected for the spurious CM motion 

_ / da [eN{a)N^{a)NM + e„{a)NA{a)NM + eMNA{a)N^{a)] , . 
^"■=° J da[<i;\'iP{a)>] ' ^ ^ 

We stress that both the nucleons and the mesons were taken into account in the evaluation 
of this projected energy. 
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IV. THE NUCLEAR CHARGE ROOT-MEAN-SQUARE RADIUS 

We now turn to the evaluation of the nuclear root-mean-square (RMS) radius in the 
formalism. Most of the measurements refer to the proton charge RMS radius so we restrict 
ourselves to that case (measurements for the neutron RMS radius are under way and are 



receiving an increasing interest ( I4])). On the other hand, since the mesons in the model 
are all neutral we have to consider just the nucleon (proton) contribution. Finally, in the 
discussion below, we consider point-particle nucleons, though nucleon form factors can be 
included without major difficulties. 

The (translationally invariant) nuclear radius operator is 

A 

R^i = J2ei{fi-Rcu? , (49) 

1=1 

where -Rcm is the center of mass coordinate and Cj is the charge of the i-th particle. The 
above operator can be rewritten, for N = Z, as: 

i=l i<j 

As for the energy calculation, the above radius operator commutes with the total linear 
momentum, but our model wave function I-?/) > is not a total momentum eigenfunction, so 
the mean radius is then given by: 

1 f da< iJ\I4,exp'^-■'\^p > 
^ Ze J da <^P\ V^(a) > ' ^ ^ 

Noting that the radius operator contains an one-body and a two-body term, the numer- 
ator of the above equation can be worked out with the help of equations like (j^6|) and (H7|) 
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respectively. 

V. NUMERICAL APPLICATIONS FOR N = Z CLOSED SHELL NUCLEI 

In order to perform applications to specific nuclei, we must solve first the a — uj model 
above described in the Hartree approximation, disregarding the CM motion effects. This is 
totally equivalent to solve the model trea ting the mesons as classical fields ||6i]. We choose to 



follow the method described in reference 13|, where both the nucleon Dirac spinors and the 
fields are expanded in three-dimensional harmonic oscillator functions, Rki{r), and treat the 
expansion coefficients as variational parameters. As we are dealing here with closed shell 
nuclei only, we have: 

N 

9niAr) = Y.^t''^^ki{r), (52) 

fe=0 

n' 

Mr) = J2ct'^RM{r), (53) 

with g and / being the up and lower radial components for the single-particle wave function. 
For the meson fields: 

Ng 

Bir) = J^CiR,oir), (54) 

k=0 

where B stands for 0o or ljo and Rki for the radial harmonic oscillator function. Those 
expansions can be introduced in the Dirac and Klein-Gordon equations and solved self- 
consistently for the expansions coefficients Ck, Ck and . After that, it is straightforward 
to implement the calculation of the energy and RMS radius as presented in the above 
sections, including the CM motion correction due to the nucleons and mesons. 

In Table I we show our results for the energy and for the root-mean-square charge radius 
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without and with the CM projection (the set of parameters for the nucleon and meson masses 



and for the couphng constants are taken from reference [ij], but disregarding the p meson 
and the eletromagnetic field). In Table II, we show the effect of the CM correction over the 
total energy, without the meson contributions, i.e., only the nucleonic degrees of freedom are 



taken in to account 3|], together with the usual harmonic oscillator approximation 13j, and 
also including the correction computed from < P\/2AM > . From Table II it is clear that 
the last two corrections are similar and not very different from the Peierls-Yoccoz correction 
without the meson degrees of freedom. Let us remember that the Peierls-Yoccoz method 
gives us not only the energy correction but also a translationally invariant wave function for 
the system. 

It is worthwhile to note that the inclusion of the mesonic contribution makes the system 
too much bounded in comparison with the case where just the fermionic contribution is 
explicitly taken in to account. However, with a slight modification of model parameters, we 
are able to obtain reasonable results for the energy and charge radius, as shown in Table III, 
in which the experimental results are also displayed. For comparison within our calculation, 
we have extracted the proton form factor contribution from the experimental charge radius 
using the prescription given in equation (6.2), ref. 13|. 

We must stress that the results shown in Table III are not obtained from a careful fitting 
of the model parameters, which should be done only after the inclusion of other mesons, as 
well as non-linear terms in the original Lagrangian. Formally, these terms can be readily 
included but then the calculations become more involved. 
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VI. CONCLUSIONS 

We have computed the center-of-mass correction in the binding energy and charge radius 
for spherical N = Z nuclei using the well known Peierls-Yoccoz projection method applied 
to the Hartree solution of the Walecka a — uj model. Although no explicit reference has to 
be made to the mesonic states in the Hartree approximation, we have chosen coherent states 
to describe meson degrees of freedom. Those states are then completely determined in this 
approximation and this allows us to obtain the nucleonic as well as the mesonic center-of- 
mass motion correction. The numerical results show a very important contribution from the 
mesons to the final binding energies and a modest but still noticeable contribution to the 
charge radius, as compared to the case where only the nucleonic CM correction is taken into 
account or to the situation where no correction is done. It is known that the Peierls-Yoccoz 
projection suffers from the so-called mass parameter problem which can be circumvented by 



using the Peierls-Thouless or the so-called variation-after-projection method 10|. Both are 



technically difficult to implement but the latter might be feasible in systems of nucleons and 



mesons, at least approximately. However, it was shown in Ref. [15|] that some observables , 
calculated in p = states, do not suffer from the Peierls-Yoccoz mass problem. We intend 
to perform a partial variation-after-projection in a restricted meson space but do not expect 
large discrepancies for oxygen and calcium, whose number of particles is already large, so 
that quantum fluctuations are expected to be smaller. 

Another important feature of our result is the fact that the CM correction, including 
the mesons, makes the system too much bounded. This is expected as long as the model 

16 



parameters were chosen to reproduce some aspects of finite nuclei witliout that correction. 
We have also shown that a few percent change in the coupling constants can bring the total 
energy and rms radius close to the experimental values, at least for the ^^O and ^°Ca cases. 
For ^He the results are still too far from the desirable using our proposed values, but this is 
true even when no CM correction is included and using the original parametrization for that 
nucleus. Furthermore, the energy correction in the '^He case is relatively large irrespective 
to the approximation used to extract the CM motion, so we believe that for this light mass 
region the projection after variation procedure may not be applicable. 

In short, we may say that, if we want to take into account the mesons in the center- 
of-mass correction applied to a relativistic model for the nucleus, a new set of parameters 
must be found in order to reproduce some basic nuclear properties as the binding energy 
and radius. Once this is achieved, it would be interesting to obtain other important nuclear 
properties as, e.g., the eletromagnetic form factors and spectroscopic factors. The model 
and the techniques explored in this paper would provide a good opportunity to obtain those 
observables. This work is in progress. 

This work was partially supported by CNPq (Brazil) and GRICES (Portugal). 



[1] K.W. Schmid and F. Griimer, Z.Phys. A337 (1990) 267. 
[2] R.R. Rodriguez-Guzman and K.W. Schmid, arXiv:nucl-th/0503059v l (2005). 
[3] P. Alberto, S.S. Avancini, M. Fiolhais and J.R. Marinelli, Brazilian Journal of Physics, 36 
(2006) 1375-1378 

[4] J.M. Eisenberg and W.Greiner, Nuclear Theory Vol.3, North Holland Publishing Company, 
1972. 



17 



[5] P. Alberto, S.S. Avancini and M. Fiolhais, International Journal of Modern Physics E14 (2005) 
1171. 

[6] B. Serot and J.D. Wallecka , Advances in Nuclear Physics 16 (1986) 1. 

[7] T.D. Lee, Particle Physics and introduction to field theory, Harwood Academic Publishers, 

New York, 1981. 
[8] M.C. Birse, Prog. Part. Nucl. Phys, 25 (1990) 1. 

[9] Shan-Gui Zhou, Jie Meng, and P. Ring, Nucl. Phys. A 694 (2001)249. 

[10] P. Ring and P. Schuck, The Nuclear Many-Body Problem, Springer, Heidelberg (1992). 
[11] J. P. Blaizot and G. Ripka, Quantum Theory of Finite Systems, The MIT Press, Cambridge, 
Massachussets, 1986. 

[12] C.J. Horowitz, S.J. Pollock, P.A. Souder and R. Michaels, Phys. Rev. C63 (2001) 025501. 

[13] Y.K. Gambhir, P. Ring and A. Thimet, Annals of Physics 198 (1990) 132. 

[14] C.J. Horowitz and B.D. Serot, Nucl. Phys. A368 (1981) 503. 

[15] E. G. Lubeck, E. M. Henley and L. Wilets, Phys. Rev. D35 (1987) 2809. 



18 



TABLE I: Ground-state energy without {E) and with (E'proj) the CM correction for the three 
double-closed shell nuclei considered in this work and the root-mean-square charge radius 
without(< >) and with (< >nio\) the same corrections. 



Nucleus 


^[MeV] 


SprojiMeV] 


< >^/^ [fm] 




^He 

160 

40Ca 


-4.85 
-94.63 
-33L32 


-68.95 
-190.67 
-420.17 


2.06 
2.59 
3.33 


1.84 
2.51 
3.28 


TABLE IL Ground-state energy, E, without the CM correction for the three double-closed shell nu- 
clei considered in this work and with the CM correction -Eproj but not considering the meson degrees 
of freedom. Also shown is the energy with the CM correction, -Eharnu calculated in the harmonic 
oscillator approximation and the energy corrected just by the subtraction of < P'^/2AAI >. 


Nucleus 


E[MeV] 


^proj[MeV] 


^harm[MeV] 


^Pl/2AM[MeV] 


^He 

160 

40Ca 


-4.85 
-94.63 
-33L32 


-18.07 
-107.87 
-342.56 


-24.22 
-106.83 
-340.31 


-16.35 
-104.92 
-339.84 



TABLE III: Ground-state energy -Eproj for the three double-closed shell nuclei considered in this 
work and charge radius < >proj '^ith the CM corrections included, compared to the experimental 
results. The figures were obtained using the values gs = 10.45 , = 13.82 and nis = 522 MeV, as 
compared to the values gs = 10.47, g^ = 13.80 and rus = 522 MeV from [l^ . 

Nucleus Sproj[MeV] ^cxp[MeV] <r2>^/J[fm] < >lil [im] 

^He -53.50 -28.30 2.01 1.57 

^^O -158.50 -127.68 2.60 2.61 

"^^Ca -339.14 -338.00 3.36 3.39 
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